6 Search Results

The structure and dynamics of confined suspensions of particles of arbitrary shape are of interest in multiple disciplines from biology to engineering. Theoretical studies are often limited by the complexity of long-range particle-particle and particle-wall forces, including many-body fluctuating hydrodynamic interactions. Here, we report a computational study on the diffusion of spherical and cylindrical particles confined in a spherical cavity. We rely on an immersed-boundary general geometry Ewald-like method to capture lubrication and long-range hydrodynamics and include appropriate non-slip conditions at the confining walls. A Chebyshev polynomial approximation is used to satisfy the fluctuation-dissipation theorem for the Brownian suspension. Wemore » explore how lubrication, long-range hydrodynamics, particle volume fraction, and shape affect the equilibrium structure and the diffusion of the particles. It is found that once the particle volume fraction is greater than 10%, the particles start to form layered aggregates that greatly influence particle dynamics. Hydrodynamic interactions strongly influence the particle diffusion by inducing spatially dependent short-time diffusion coefficients, stronger wall effects on the particle diffusion toward the walls, and a sub-diffusive regime-caused by crowding-in the long-time particle mobility. The level of asymmetry of the cylindrical particles considered here is enough to induce an orientational order in the layered structure, decreasing the diffusion rate and facilitating a transition to the crowded mobility regime at low particle concentrations. Our results offer fundamental insights into the diffusion and distribution of globular and fibrillar proteins inside cells.« less

Here, we report a computational strategy to obtain the charges of individual dielectric particles from experimental observation of their interactions as a function of time. This strategy uses evolutionary optimization to minimize the difference between trajectories extracted from the experiment and simulated trajectories based on many-particle force fields. The force fields include both Coulombic interactions and dielectric polarization effects that arise due to particle-particle charge mismatch and particle-environment dielectric contrast. The strategy was applied to systems of free falling charged granular particles in a vacuum, where electrostatic interactions are the only driving forces that influence the particles' motion. We showmore » that when the particles' initial positions and velocities are known, the optimizer requires only an initial and final particle configuration of a short trajectory in order to accurately infer the particles' charges; when the initial velocities are unknown and only the initial positions are given, the optimizer can learn from multiple frames along the trajectory to determine the particles' initial velocities and charges. While the results presented here offer a proof-of-concept demonstration of the proposed ideas, the proposed strategy could be extended to more complex systems of electrostatically charged granular matter.« less

Molecular simulation has emerged as an essential tool for modern-day research, but obtaining proper results and making reliable conclusions from simulations requires adequate sampling of the system under consideration. To this end, a variety of methods exist in the literature that can enhance sampling considerably, and increasingly sophisticated, effective algorithms continue to be developed at a rapid pace. Implementation of these techniques, however, can be challenging for experts and non-experts alike. There is a clear need for software that provides rapid, reliable, and easy access to a wide range of advanced sampling methods, and that facilitates implementation of new techniquesmore » as they emerge. Here we present SSAGES, a publicly available Software Suite for Advanced General Ensemble Simulations designed to interface with multiple widely used molecular dynamics simulations packages. SSAGES allows facile application of a variety of enhanced sampling techniques—including adaptive biasing force, string methods, and forward flux sampling—that extract meaningful free energy and transition path data from all-atom and coarse grained simulations. A noteworthy feature of SSAGES is a user-friendly framework that facilitates further development and implementation of new methods and collective variables. In this study, the use of SSAGES is illustrated in the context of simple representative applications involving distinct methods and different collective variables that are available in the current release of the suite.« less

Here, we studied the transport of room-temperature ionic liquids (RTILs) through charged conical nanopores using a Landau-Ginzburg-type continuum model that takes steric effect and strong ion-ion correlations into account. When the surface charge is uniform on the pore wall, weak current rectification is observed. When the charge density near the pore base is removed, the ionic current is greatly suppressed under negative bias voltage while nearly unchanged under positive bias voltage, thereby leading to enhanced current rectification. These predictions agree qualitatively with prior experimental observations, and we elucidated them by analyzing the different components of the ionic current and themore » structural changes of electrical double layers (EDLs) at the pore tip under different bias voltages and surface charge patterns. These analyses reveal that the different modifications of the EDL structure near the pore tip by the positive and negative bias voltages cause the current rectification and the observed dependence on the distribution of surface charge on the pore wall. The fact that the current rectification phenomena are captured qualitatively by the simple model originally developed for describing EDLs at equilibrium conditions suggests that this model may be promising for understanding the ionic transport under nonequilibrium conditions when the EDL structure is strongly perturbed by external fields.« less

Large classes of materials systems in physics and engineering are governed by magnetic and electrostatic interactions. Continuum or mesoscale descriptions of such systems can be cast in terms of integral equations, whose direct computational evaluation requires O(N²) operations, where N is the number of unknowns. Such a scaling, which arises from the many-body nature of the relevant Green's function, has precluded wide-spread adoption of integral methods for solution of large-scale scientific and engineering problems. In this work, a parallel computational approach is presented that relies on using scalable open source libraries and utilizes a kernel-independent Fast Multipole Method (FMM) tomore » evaluate the integrals in O(N) operations, with O(N) memory cost, thereby substantially improving the scalability and efficiency of computational integral methods. We demonstrate the accuracy, efficiency, and scalability of our approach in the context of two examples. In the first, we solve a boundary value problem for a ferroelectric/ferromagnetic volume in free space. In the second, we solve an electrostatic problem involving polarizable dielectric bodies in an unbounded dielectric medium. Lastly, the results from these test cases show that our proposed parallel approach, which is built on a kernel-independent FMM, can enable highly efficient and accurate simulations and allow for considerable flexibility in a broad range of applications.« less